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Abstract 

We consider an example of a quantum algorithm from the point of view 
of the de Broglie-Bohm formulation of quantum mechanics. For concrete- 
ness we look at two particular implementations: one using spin-^ particles 
as described by a simple model due to Bell and the other as energy states 
in an infinite-well potential. We extend the analysis to encompass a com- 
plete set of quantum gates. We conclude by discussing the relevance of our 
investigations for the debate about the origin of the quantum-computational 
speed-up. 

1 Introduction 

Exploiting quantum phenomena for performing computation allows an exponen- 
tial speed-up in particular types of computational tasks, such as the factorization of 
large numbers (Shor, 1997, [[TJ). While the field of quantum computation has seen 
dramatic advances in the last two decades, no consensus has been reached regard- 
ing the question of exactly which aspects of quantum physics allow the speed-up 
relative to computations performed by classical means. Numerous proposals have 
been made, including entanglement (Josza, in Huggett, 1998, [|2)), the wave-like 
behaviour of quantum systems (Mermin, 2007, 0) and the existence of "many 
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worlds" or "parallel universes" (Deutsch, 1985, 1997 [|4j[5j and in Saunders et al., 
2010, §|). 

We consider it instructive to look at quantum computing from the perspective 
of another realist theory, de Broglie-Bohm Pilot- Wave Theory (PWT). Whether 
or not Many-Worlds Theory (MWT) has the explanatory power that some of its 
proponents say it has with regards to quantum computing, and whether or not it 
is the only "interpretation" of quantum mechanics with that power, is a debate 
that we will not engage in here. Either way, a fair assessment of the issue relies 
on a full understanding of quantum computation from the point of view of any 
contending theory. Remarkably, while PWT has received growing interest during 
the last twenty years or so, the present author has been unable to find an analysis 
of even a simple quantum algorithm in the context of the theory. 

We hope to remedy this situation here by examining the simplest non-trivial 
example of a quantum algorithm, namely the so-called Deutsch algorithm[j] In 
order to be able to compute de Broglie-Bohm trajectories we will have to make a 
choice concerning a particular physical model for the qubits. In this paper we ex- 
amine two such models, spin-^ particles as described by Bell's simple spin model 
(Bell, 1987, [7 1), followed by a representation of |0) and |1) as the two lowest en- 
ergy levels of an infinite square well. We then extend this analysis to a complete 
set of gates, i.e. a set of gates that is sufficient to implement any unitary evolution 
for an arbitrary number of qubits. Thus the results of this paper show explicitly 
how PWT can describe elementary quantum computations. Some possible foun- 
dational implications are sketched in the conclusion. 

2 Review of the de Broglie-Bohm formulation 

The de Broglie-Bohm pilot-wave formulation of quantum mechanics is a non- 
local, deterministic hidden-variable theory. Here we will focus on the non-relativistic 
case, since this will be relevant for our analysis of the quantum algorithms below. 
Its ontology contains a wavefunction ip, evolving in configuration space accord- 
ing to the Schrodinger equation in the position basis as familiar from the standard 

'This is the single-qubit case of the more general "Deutsch-Josza" algorithm. 
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Copenhagen interpretation, 



i—il>(x,t) = Hip(x,t) 



(using h — 1), and an actual configuration (the "particle" or "corpuscle") forming 
a trajectory in configuration space according to the "guidance" equation 

d _ j(x,t) 
-x - 



dt \^(x,t)\ 2 

where j (x, t) is the standard probability current, implicitly defined by 

||^(f,t)| 2 + V-j(f,t) = 0. 
For an n-particle system with a Hamiltonian 

n 

8=0 l 

where m 8 denotes the mass of the 2th particle and Vj the derivative with respect 
to its spatial degrees of freedom, the guidance equation may be simplified to 

^-Si = —VS(x,t) = —3m(Vihiib(x,t)) 
dt rrii rrii 

with Xi denoting the particle's position in the ith particle's subspace of the full 
configuration space. S(x, t) is the phase of ip(x, t). 

A measurement in PWT is the interaction between the system and a measure- 
ment apparatus according to an effective Hamiltonian^] 

d 

H m eas l(lQ — , 

dy 

where y is a variable for the degree of freedom of the "pointer" or "dial" on the 
apparatus from which we may read off the "outcome" of the measurement, Q 
is a unitary operator associated with the particular measurement and agrees with 
the corresponding measurement operator in the standard formulation, and a is an 

2 The actual Hamiltonian would likely be considerably more complicated given that most ap- 
paratuses are made up of very large numbers of atoms. 
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appropriate coupling constant. The actually observed outcome then depends on 
the position of the configuration- space particle (see Holland, 1993, (5J). 

Instead of using the guidance equation we may choose to adopt a second- 
order equation analogous to the Hamilton- Jacobi equation with an extra term (see 
Bohm, 1952, [(9)), 

i=l 

where 

1 V?|#?,f)| 



) _ | Inii \ip(x,t)\ 

is known as the "quantum potential" due to its role in the equation being analogous 
to the potential V. Note however that U, unlike V, is calculated from the wave- 
function and strongly time-dependent. In this paper, we will use the first-order 
formulation of PWT. 

If many copies of the same system are available, then this ensemble may be de- 
scribed by density function p(x, t) in configuration space. If p(x, t) = \ijj(x, t)\ 2 , 
then we say that the ensemble is in "quantum equilibrium" (Valentini, 1991a, 



[10 1). It is straightforward to show that the observed frequency of outcomes 
matches the predictions of the Born Rule of standard quantum mechanics, which 
states that the possible outcomes of a projective measurement associated with 
operator Q are the operator's eigenvalues and the probability of obtaining the ith 
eigenvalue A, given the system's initial state \ip) is | ( ip\ i) | 2 where \i) is the eigen- 
state of Q associated with eigenvalue Aj. 

If p(x, t) 7^ \i)(x, t)\ 2 , then the ensemble is in non-equilibrium and outcomes 
of measurement show statistical deviation from the Born Rule. For ensembles 
with certain non-equilibrium density functions faster-than-light signalling and vi- 
olations of the uncertainty principle are possible (Valentini, 1991b, (10J). How- 
ever, subject to certain conditions the ensemble will approach equilibrium as it 
evolves and a theorem analogous to the classical coarse-graining H-theorem may 
be derived, explaining the empirical adequacy of the Born Rule. This paper is pri- 
marily concerned with the equilibrium case, although much of our analysis applies 
to both equilibrium and non-equilibrium. 
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3 Deutsch's algorithm and generating Hamiltonians 



The simplest known non-trivial quantum algorithm is the so-called Deutsch algo- 
rithm. We are given a black box (the "oracle") that implements one of the four 
possible functions / : {0, 1} — > {0, 1}. Which of these it implements is unknown 
to us. Enumerate the possible functions as follows: 



/o(0) 


= 0; 


/o(l) 


= 


/i(0) 


= 0; 




= 1 


/ 2 (o) 


= i; 




= 


/a(0) 


= i; 


m) 


= 1 



We are specifically tasked to determine if / is constant (i.e. / G {fo, /s}) or 
balanced (i.e. / G {/i, /2})F] Classically this task requires two uses of the ora- 
cle, once with input '0' and once with '1' since neither by itself narrows down 
the set of possible functions sufficiently. Exploiting quantum mechanics (exactly 
what aspects of quantum mechanics those are is the debate to which we hope to 
contribute through this analysis) however only one use of the oracle is required. 

To see how, we first make the mechanics of the oracle more explicit: Two 
qubits are passed through the oracle, the first of which serves as providing the 
argument of the function / implemented by the oracle (call this the "data qubit", 
following the terminology of Nielsen & Chuang (2000, p~TJ)) and it does not 
change state as it passes through. The second ("auxiliary") qubit carries the out- 
put: At the end of the oracle interaction it will be in a state \q © f(p)), where 
q,p G {0, 1} denote the input states of the auxiliary and data qubits respectively 
and '©' is used to denote addition modulo 2 . The action on superposition states 
follows from linearity. 

To achieve the stated task, begin with the qubits in the combined state |0) d |1) , 
where the subscripts d and a identify the data and auxiliary qubits respectively. We 
will express all two-qubit states in the product basis Sj d ^S) a , i.e. the tensor product 
of the two two-dimensional single-qubit Hilbert spaces with basis { |0) , 1 1) }. First 

3 For other tasks, such as to determine / fully or to determine whether or not / = /i, say, 
quantum mechanics is not known to allow any computational speed-up. 
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apply a Hadamard gate H 



1 

V2 



1 



to each qubit, followed by the action 



of the oracle, which may be stated in compact 2x2 form as 

Uf = 



<W(o)l + di,f(o)X 

+ Sij(i)X 



(1) 



1 

1 



denotes the first Pauli 



where S it j is the Kronecker-delta symbol and X - 

matrix. Finally act with another Hadamard gate on the data qubit and perform a 
measurement on it in the {|0) , |l)}-basis. It is easy to verify that an outcome of 
implies that / is constant and an outcome of 1 that / is balanced. For a more 
detailed description we refer the reader to Nielsen & Chuang (2000, fTTJ). 

As preparation for our analysis we need to find the generating Hamiltonians 
for the unitary evolutions of the quantum gates that constitute the circuit, namely 
the Hadamard gates and the two-qubit evolution Uf. The Hadamard gate is gen- 
erated by a Hamiltonian of the form 



H 



Had 



T, 



7T 



Had 



where T Had denotes the duration of the action of H Had and 1 is the identity oper- 
ator, as may be readily verified. Thus tJ— may be understood as the parameter 

Had 

determining the "strength" of the acting Hamiltonian. 

The block-diagonal form of eq. [T] makes the task of finding the generating 
Hamiltonian for the oracle action straightforward. It is: 



H f T, 



f 1 or 



$L,/(0) ■ \{ X 





) o 

<W f(*-i) 

where analogously to above T or . denotes the duration of the oracle's action. 



(2) 



4 Bell's spin toy model 

Let us first analyse the Deutsch algorithm as it might be implemented by a toy 
model due to Bell (1987, 0). In the context of PWT the model was first stud- 



ied by Valentini (2002, [ 12 1). The model represents such a particle as a complex 
two-vector without spatial dependence. The only spatial degree of freedom is that 
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of the pointer of the spin measurement apparatus. Note that spin is not a "be- 
able", that is, there are no degrees of freedom of the configuration space particle 
associated with spin. 

Deutsch's algorithm uses two qubits and hence we need two spin-| particles, 
in addition to a single measurement device, characterised by a variable y, that 
is required for the final step of the algorithm. Initially the wavefunction of the 
measurement device is a narrow wave packed 4>o(y) centred at y = 0. The total 
initial wavefunction is then given by 

1pmn(y,to) = <f>o(y)d m (to)a n (t ), 

where we define d m (t) and a n (t) to be the wavefunctions (i.e. two-vectors) of the 
data and auxiliary qubit respectively and t is the point in time when the compu- 
tation is set into motion. 

Calculating the evolution of the pilot wave is then straightforward. Given 
that the wavefunctions of the two particles are completely specified by a two- 
component vector of unit length, the two states of the abstract description provided 
in the previous section can be directly identified with the two eigenstates of Z = 
1 

o -l J : 

\0} d ^d m = ( J ) , |l) d *+d m = ( J ) , 

\ / m V / m 

|0} a ^O„= ( J ) > \ 1 )a^ a n= ( ^ ) • 

The evolution induced by the two Hadamard gates is given by 

4>(y)d m (to)a n (t ) ->■ H mp H nq (J)o(y)dp(t )a q (t ) 

and for input states d m = ( J J and a n = ( ^ J this corresponds to 



/ m \ / n ' \ / m \ / n 



which we choose to rewrite using pair labels in a single four-component vector, 
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i.e. as 



>(</) 



with m, n G {0, 1}. 

The action of the oracle is then 



/ 1 \ 

-1 

1 

V- 1 / 



(mn) 



My) 



( i \ 

-i 

i 

V - 1 J 



-> 2^{y) U f (mn)(pq) 



(mn) 



( 1 \ 

-1 

1 

V- 1 / 



(pq) 



Denote the time when the final measurement occurs by t meas . The pilot wave 
at a time t meas — e just prior to the measurement is then given by 



j ' i 



Ipmn (?/ j tf. 



My) 



My) 



V2 



V2 



( 1 \ 

-1 



V J 

( \ 


1 

V - 1 J 



if/(0) = /(l) 



if/(0)^/(l). 



We note that up to this point the qubits have not interacted with the measurement 
apparatus, whose wave function 4>(y) is still unchanged. Hence the ensemble 
equilibrium distribution is unchanged too. 

Following standard de Broglie-Bohm measurement theory (Holland, 1993, 
(8), Bell, 1987, [7|, Bohm, 1952, (5)), the interaction Hamiltonian for the sub- 
sequent measurement in this model is given by 

, d 



H 



Bell 



-ig(Z 



dy 



yielding for a time immediately after the measurement 



j L meas + e) = ^mn{y ~ 2#e, 
i>mn(y,t meas + e) = ^mn{y + 2ge,t 

r, 



-e) 



if/(0) = /(l) 
if/(0)^/(l). 
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Thus, if the initial wave packet of the pointer is sufficiently narrow, we can un- 
ambiguously read off whether / is constant or balanced. We note that this result 
does not strictly rely on being in quantum equilibrium (p = |t/>| 2 ), but any ensem- 
ble distribution would suffice that evolves in such a way that the two packets of 
non-zero density separate with the separation of the pilot- wave packets. 

If we wish we can also trace the evolution of the pilot wave during the times 
"inside" a gate. For example, for the Hadamard gate we found the Hamiltonian 



Setting t — when the gate is "switched on", we may consider the state of the 
wave at times t — aT where a e [0,1]. The evolution operator is 



For a = 1 this reduces to the complete Hadamard gate and for a = to the 
identity. The wavefunction of the pointer (and hence the ensemble distribution) 
is naturally unchanged while passing through the gate in this model for spin. As 
such, investigating the evolution "within" a gate is only of marginal interest here, 
although it plays a greater role in our treatment of the infinite- well model.. 

We have discussed the evolution of the ensemble (unchanged until the mea- 
surement, then moving with the shift of the wave packet of the pointer) but have 
not yet calculated any actual trajectory. Prior to the final measurement the wave- 
function of the pointer and that of the qubits are entirely independent. Hence for 
this period (t < t meas ) we can write the total wave function as 



where r(t) is the normalised four-component state vector of the two qubits and 
4>(y, t) is the wavefunction of the pointer, such that <f)(y, t) = 4>o(y) for t < t meas . 
The Hamiltonian H Gate of any of the gates acts only on r(t) and is expressible as 
a Hermitian matrix. Hence for the current we find unsurprisingly 



and so if we continue to neglect the free evolution of the pilot-wave, the equation 
of motion for t < t rneas is simply 



H Had T=-[H-l]. 




il>{y,t) = 4>{y,t)r(t) 



j(y,t) 



d_ 

di 



My,t)\ 2 = o 



dy 
dt 



= 0. 
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When the measurement is carried out, the equation becomes 

dy , 
dt =±9 - 

So not only the wave packet and density distribution as a whole moves with 
velocity g in the positive or negative y-direction depending on the oracle function 
/, but in fact the configuration particle itself moves with this velocity. The trajec- 
tories of the ensemble are exactly parallel. The direction of the trajectories depend 
on the wavefunction. From this we may wish to conclude that most of the "work" 
of the computation is done by the wavefunction and the actual trajectory has little 
to do with it. 

5 An infinite well model 

We will now consider a model in which qubits are realised by the ground state 
and first excited state of particles in an infinite well, colloquially also known as 
particles in a box. In particular, to keep the description as simple as possible we 
will represent qubits by the states of a 1 -dimensional well. Since the Deutsch 
algorithm requires two qubits, we will either need two such particle-well systems, 
or a single 2-dimensional one. 

Assume that the box has a length L — 1 and is described by coordinates 
x G [0, 1]. The ground and first excited state in the position basis are then given 
by 

(x\ 0) = (x\ ipi) = V2 sin(7rx), (x\ 1) = (x\ ip 2 } = V% sin(27nr). 

The discrepancy of index labelling between the wavefunction (^1,^2) and the 
state kets (|0) , 1 1)) arises from our desire not to deviate from conventions used in 
the general literature. A general state \ip) = a \-ipx) + b \ip 2 ) evolves into 

\^(t)) = e- l(£l+S2) * /2 (ae-^|Vi) + &e iarf |^ 2 )) 

with 00 = \{Ei — E 2 ), which in abstract vector notation corresponds to a state 
( °. ) evolving according to a Hamiltonian Hf ree = uZ and so its time evolution 
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IS 



a 
b 



—> R z {2ui) ( ° ^, where i? z (#) = e iez / 2 , plus an overall phase 

rotation. This means we can implement a Z-gate, for example, by simply waiting 
for the correct amount of time cut = tt. 

Recall that any single-qubit gate can be written in the form 

U Gate = e ia R z (/3)RAl)Rz(S), 

where R x (9) and R z (x) correspond to rotations in the Bloch sphere by angles 9 
and x around the x and z-axis respectively and that geometrically any arbitrary 
rotation may be implemented by three sequential rotations round these axes. e ia 
fixes the overall phase, which in our example may be implemented using the time 
evolution of the overall phase factor e - l ( E i+ E 2)t/^_ 

Performing rotations R x (9) may be achieved by perturbing the potential (and 
hence the Hamiltonian) inside the well. In particular, adding the term 



8V(x) 



9vr2 
16 



1 



x — 



(see e.g. [ 1 1 ] p. 280) we find that its matrix elements are given by 

(tp n \ 5V \ip m ) -- 



1 

1 



which can be easily verified by performing the integration explicitly in the position 
basis. This Hamiltonian thus generates the required rotations R x around the x-axis 
in Bloch space. 



The Hadamard gate H = ^= 



1 1 
1 -1 



for example, is given by 



H = e^R z 



Rx 



Rz 



corresponding to a free evolution for a period At = followed by a perturbation 
of the potential by 5V for a period At = f and finally followed by another period 
of free evolution. We justify the neglect of the free evolution during the second 
time interval (the X-rotation) by choosing u to be small (e.g. by choosing a large 
particle mass m). The overall phase will be relevant in interactions with other 
qubits, but in the Deutsch algorithm there is only one such interaction (the oracle) 
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and hence the phase may be fixed by preparing the correct states at the right time 
to achieve the right phase factor for the interaction. 

We will now turn our attention to the interaction of the oracle. We will con- 
sider specific choices /. / is trivial to implement (no evolution at all) and f 3 is 
simply a NOT-gate in the auxiliary qubit, which we can implement by perturbing 
the potential as above for an interval At = |. In neither of these cases is any inter- 
actions between the qubits required. However, for the cases of balanced functions 
fi, i.e. for fi and f 2 , the situation is different and an actual two-qubit gate needs 
to be constructed. 

We will only consider f 2 here. A construction of fi is easily obtainable by 
modifying the description here, or by using the f 2 evolution exactly, but placing 
an additional NOT-gate before and after the oracle in the data qubit. We implement 
the oracle by perturbing the potential by an additional term U (x, y), where x and 
y denote the position coordinates of the data and auxiliary qubits respectively. 
We remark that if physically we have a single particle in two dimensions, such a 
perturbation may be easily visualised. For two particles in one-dimensional wells, 
this term corresponds to a (possibly quite complicated) coupling between them. 
Equation [2] for the case f 2 reads 



- 1 or- 1 - 1 J2 %J 



i(X-l) 




(3) 



Choosing for simplicity the time interval during which the gate acts to be 
T or = §. U(x, y) must have the following matrix elements: 

(Tpk\ d (tpm\ a U(x,y)\'ip n ) a \'il; l ) d = (X-l) mn fork = 1 = 1 
(^k\ d (^m\ a U(x,y)\i) n ) a \i/;i) d = otherwise 

with k.l.m, n E {1, 2}. This may be achieved with a potential U(x, y, ) inside the 
well of the form 

9n 2 , 1. 



U (x, y) = (A + B cos x + Cx cos x) 



(y-s)-i 



16 vw 2 



where A, B and C are appropriately chosen constants. Appropriate numerical val- 
ues are A = p, B = —^n 2 and C = |f§7r 2 , although other linear combinations 
of functions of x are also possible. The matrix elements {ip P \ d (ip m \ a U(x, y) \i>n) a \4>q) c 
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then satisfy the requirements given above as may once again be verified by explicit 
calculation in the position basis. 

The final measurement is an interaction between the data qubit and the pointer 
described by a position variable z and whose wavefunction starts as an initial 
narrow wave packet centered at z = 0. An appropriate interaction Hamiltonian is 
given by 

Hmeas = d(-id x ) 2 (-»£?*) , 

where a is a coupling constant. From PWT measurement theory we know that as 
a result the ensemble density shifts by a5t ■ n 2 ii 2 for n = 1,2 corresponding to 
the two energy states, and 5t denotes the duration of the interaction between qubit 
and apparatus. The projection of the trajectories onto the configuration space 
dimension of the pointer are all exactly parallel. 

We see that we can fully account for the outcome of the computation by only 
considering the wavefunction up to the point of measurement and then appeal to 
standard PWT measurement theory. As such, we have also automatically analysed 
the evolution of the ensemble, at least in quantum equilibrium. A corresponding 
analysis for quantum non-equilibrium would possibly present interesting future 
work, although we can already conclude that if -0 = — > p = (i.e. the density 
function only has "support" where the pilot-wave is non-zero)j^]then the computa- 
tion works exactly the same nevertheless: The measurement outcome is determin- 
istic in every sense (as follows from the quantum formalism) and only depends on 
the choice of oracle function /, hence any trajectory whose corpuscle is at some 
time inside the configuration space area where ip ^ corresponds to the correct 
computational output. This conclusion does not apply to other, non-deterministic 
algorithms. 

The only task left before us is to calculate trajectories of the configuration 
space particle (in the x-y subspace of configuration space) during the action of 
the gates. Note that the aspect most crucial for the trajectory is its projection 
onto the dimension corresponding to the degree of freedom of the pointer of the 

4 For this not to be the case the ensemble density would have to be non-zero somewhere outside 
the box. The fact that the phase is undefined in such areas will likely be problematic for a coherent 
PWT description for such "extreme" non-equilibrium. 
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measurement apparatus, which we have already discussed. We may wish to in- 
vestigate the trajectories of the qubits themselves for completeness. Since we 
are dealing with systems of spinless particles whose Hamiltonian is of the form 
H W eii = — + V, the guidance equation is the usual mx = VS. We are par- 
ticularly interested in how the trajectories of the particles differ between the cases 
of a constant and a balanced oracle function /. We recall that the data qubit enters 
the oracle in a state \+) d = sin(7rx) + sin(27rx)). If / is constant it will return to 
that state when the oracle action stops and if / is balanced its final state will be 
|— ) d = sin(7r:r) — sin(27rx)). The corresponding quantum equilibrium ensemble 
density functions in the ^-dimension are depicted in figure [T] below. 




x position x-position 



Figure 1: Ensemble densities p(x) = \ (x\ +) | 2 = | sin(7nr) + sin(27rx)| 2 and p(x) = 
| (x\ — ) | 2 = | sin(7rx) — sin(27rx)| 2 for an infinite well of length L = 1. 

This already highlights the key difference between the types of trajectories 
found in the two cases: In the constant case, the particle ensemble returns to its 
initial distribution of x-positions, in the balanced case it does not. For both a 
constant and a balanced oracle function, the auxiliary qubit evolves from |— ) a to 
|— ) a (i.e. its initial and final states are the same) and so the ensemble returns to its 
initial y-distribution. 

We will now take a glimpse at the trajectories themselves. Since, for pi- 
lot waves in non-trivial superpositions of the two eigenfunctions sin(7rx) and 
sin(27rx), the phase S is given by the inverse tangens of a complicated function 
of x, y and t and the trajectories cannot be found analytically, we analyse these 
numerically. Here we will only consider the evolution of the configuration particle 
during the action of the oracle rather than the entire algorithm. We will attempt 
to compare the two specific cases f and f 2 , i.e. one constant and one balanced 
oracle function /. 
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First, however, consider the free evolution of the qubits in the absence of any 
gates. The phase S of the states |0) and |1) is just given by —%E\t and —iE 2 t 
respectively and hence VS = in both cases. The particle is at rest. This is not 
the case for superpositions. The phase of some general state 

|^) = ae~ iuit V2sm(7rx) + be iu>t V2sm(2irx) 

(up to some additional overall phase factor independent of x), is given by 

f — asminx) + 6sin(27rx) 1 

S = arctan < tan(wt) — r — >. 

[ asm{7ix) + bsm{2iTX) J 

Its gradient VS determines the particle trajectories and results in a one-dimensional 
oscillatory behaviour in xy-space. 

Let us now investigate the oracle. For the constant function f , the corpuscle 
remains at rest (neglecting its free evolution), from which the ensemble evolution 
anticipated in figure [T] follows trivially. Let us contrast this to the evolution gen- 
erated by the corresponding Hamiltonian, given in equation [3] One result to note 
immediately is that the particle does not move at all in the y-dimension, no mat- 
ter what its initial values for x and y are. Thus its ensemble density evolution is 
again recovered naturally. The x-evolution is more interesting. The results of the 
numerical modelling can be seen below in figure [2} The initial distribution of the 
particle positions corresponds to (coarse-grained) quantum equilibrium. 




0.2 0.4 0.6 O.B 1 1.2 1.4 1.6 1.6 

time 



Figure 2: Trajectories in the x-dimension for an oracle implementing f 2 
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The trajectories are consistent with our expectation. The ensemble density 
shifts in x-space from its initial to its final distribution in relatively simple paths. 
Recall that in the case of Bell's spin model spatial motion was only necessary for 
the pointer in order to account for the computation of the Deutsch algorithm. 

An infinite well is not a realistic realisation fo qubits and parameters for the 
numerical simulation were chosen arbitrarily. Future work should include more 
realistic realisations of the qubits where parameters such as the particle mass can 
be chosen based on physical considerations. 



6 Description of a complete set of gates 

A de Broglie-Bohm analysis using Bell's model for any quantum algorithm is 
relatively trivial: Abstract |0), |1) states map straightforwardly onto spin states 
and the only degrees of freedom for the particle are those of the measurement 
apparatus, which have already been discussed, although in general multiple qubits 
may be measured. 

Let us therefore instead extend our analysis of the infinite-well model. It has 



been shown (DiVincenzo, 1995, (T3J) that three gates are together sufficient for 

qubits. One possible such set is 

(both single-qubit gates) and the 



any desired unitary evolution for any number of qubits. One possible such set is 

"1 

e^ /4 



the Hadamard gate, the f-gate T = 

(two-qubit) controlled-NOT gatej^] 

The last of these we have already discussed in detail: The function f 2 of 
Deutsch's algorithm corresponds exactly to a controlled-NOT gate (possibly up 
to a NOT-gate before and after the controlled-NOT, depending on one's preferred 
definition). 

The | -gate is generated by a Hamiltonian 

Hk = — - -(Z-l), 

8 Tk 8 

8 

where Te is the interval during which this Hamiltonian is "switched on". An 
analysis of the f-gate is simple if the initial state is one of the basis states. In that 



5 As discussed above, any single qubit gate can be implemented using only z and a; -rotations, 
and an overall phase. Hence these together with the CNOT-gate would also form a complete set 
for constructing unitary gates with any number of qubits. 
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case the evolution is limited to the overall phase and the configuration particle 
does not move. For other initial states the evolution corresponds just to the free 
evolution. The phase factors of the two basis states rotate in opposite direction, in 
addition to an irrelevant change in the overall phase. The resulting behaviour of 
the particle is an oscillation inside the well that depends on its initial configuration. 
Numerical analysis for a small equilibrium ensemble of such particles has yielded 
the trajectories given in fig. [3} 




t 1 1 1 1 1 r 

1 2 3 4 5 6 

Time 



Figure 3: Sample trajectories for the free evolution, i.e. with an unperturbed 
infinite- well potential. 

Finally recall the Hadamard gate discussed in section [3] Consider the concrete 
evolution example |1) — > H |1) = r^(|0) — The corresponding evolution of 
the equilibrium distribution in the infinite well therefore takes the form 



\/2sin(27ra;; 



V2 



sin(7ra;) - \/2sin(27ra;) 



Fig. [4] shows the evolution of 20 trajectories constituting a coarse-grained equilib- 
rium ensemble. 
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time 



Figure 4: Sample trajectories for the Hadamard gate in the infinite-well model, 
with initial state |1) corresponding to the well's first excited state. 

Of course, here we have only displayed concrete trajectories for chosen initial 
states. Other initial states could be used, leading to different trajectories, although 
ultimately no fundamentally new insight would be gained from a more compre- 
hensive analysis in this respect. 

While we have provided a description of trajectories for a set of gates that 
together are sufficient to construct a circuit equivalent to any chosen multi-qubit 
evolution, it might nonetheless be instructive to look at another algorithm as a 
whole, bearing in mind that our purpose is not to show that PWT is able to account 
for quantum computation (that already follows from the fact that equilibrium PWT 
reproduces standard quantum mechanics empirically), but to describe how it does 
so in order to facilitate a discussion of how the relative speed-up of quantum 
computation can be explained. However, we will leave this to possible future 
work. 
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7 Conclusion and Discussion 



We have presented an explicit example of a simple quantum algorithm and a de- 
scription of a complete set of gates in the context of PWT. While the theorem 
that equilibrium PWT reproduces all measurement results of standard quantum 
mechanics exactly already guarantees that quantum computers work in PWT, it is 
nonetheless instructive to see this done explicitly. 

Quantum computation provides insightful and interesting examples of quan- 
tum evolutions and measurements. The original motivation for this paper was to 
understand how such processes could be described in PWT. Having done so, one 
may ask what implications this work may have for the interpretation of quantum 
computing and quantum mechanics generally. While in our view fundamental 
questions concerning quantum theory are not really tied to quantum computing, 
even if it provides striking examples, some authors have claimed otherwise and 
it therefore seemed worthwhile to study quantum computing from the perspective 
of PWT. 

What conclusions can we draw? We have seen that the essence of the compu- 
tation lies with the evolution of even in PWT. This was particularly evident in 
our analysis in terms of Bell's model of spin. Even so, it should be emphasised that 
the pilot-wave account is conceptually and physically distinct from that of many- 
worlds theory, despite some claims to the contrary. In particular, its account of 
measurement is radically different and outcomes are accounted for in terms of an 
actual configuration, not a global object such as the wavefunction. ip is not inter- 
preted as constituting many worlds, but simply as a complex field in configuration 
space. Everettians may refer to the often cited "many worlds in denial" argument 
(see Deutsch, 1996, [ 14], and Brown & Wallace, 2007, (13)) and argue that PWT 
involuntarily includes the structure of many worlds. However, this argument is 
contentious (e.g. Valentini, in Saunders et al., 2010, j6j) and furthermore there 
are, arguably, open issues with the claim concerning the emergence of the many- 
world structure in the Everettian picture^] Either way, it is our hope to contribute 
to this ongoing debate with the technical analysis provided in this paper. 



For a discussion of the author's views we refer the reader to Roser ( 1 16 ). 
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